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Abstract 

We investigate strong coupling effects on the superfluid phase transition in 
a gas of Fermi atoms with a Feshbach resonance. The Feshbach resonance 
describes a composite quasi-Boson, which can give rise to an additional pair- 
ing interaction between the Fermi atoms. This attractive interaction becomes 
stronger as the threshold energy (2v) of the Feshbach resonance two-particle 
bound state is lowered. In a recent paper, we showed that in the uniform Fermi 
gas, this tunable pairing interaction naturally leads to a BCS-BEC crossover 
of the Nozieres and Schmitt-Rink kind, in which the BCS-type superfluid 
phase transition continuously changes into the BEC-type as the threshold en- 
ergy is decreased. In this paper, we extend our previous work by including 
the effect of a harmonic trap potential, treated within the local density ap- 
proximation (LDA). We also give results for both weak and strong coupling to 
the Feshbach resonance. We show that the BCS-BEC crossover phenomenon 
strongly modifies the shape of the atomic density profile at the superfluid 
phase transition temperature T c , reflecting the change of the dominant parti- 
cles going from Fermi atoms to composite Bosons. In the BEC regime, these 
composite Bosons are shown to first appear well above T c . We also discuss 
the "phase diagram" above T c as a function of the tunable threshold energy 
2u. We introduce a characteristic temperature T*(2v) describing the effective 
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crossover in the normal phase from a Fermi gas of atoms to a gas of stable 
molecules. 

PACS numbers: 05.30.Jp, 03.75.Fi,67.90.+z 
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I. INTRODUCTION 



The search for the superfluid phase transition in a trapped gas of Fermi atoms is one 
of the most exciting problems in current physics. |lj This involves the formation of BCS 
Cooper-pairs made up of degenerate Fermions of two different species ( "spin up" and "spin 
down"). Degenerate Fermi gases have been realized in 40 K [0] and 6 Li. Very recently, 

a Feshbach resonance was observed in 40 K working with the two hyperfme states |9/2, —9/2) 
and |9/2, — 7/2). || In the Feshbach resonance, a quasi-molecular Boson is formed, which 
can produce an attractive pairing interaction between two Fermi atoms of opposite spins. 

As a result, one expects a BCS-type superfluidity with a high phase transition 
temperature T c due to this strong attractive interaction. |§ However, such strong interaction 
is also known to enhance fluctuations, which can strongly suppress T c predicted by the 
usual mean-field BCS theory. Thus the inclusion of these Cooper-pair channel fluctuations 
is crucial in considering the high-T c superfluidity induced by an atomic Feshbach resonance. 



In a previous paper, JTTJ we have discussed how this fluctuation contribution affects 
the superfluid phase transition temperature in a uniform Fermi gas, extending the theory 



developed by Nozieres and Schmitt-Rink: fll2] , |13],|14|,|15[ We assume that iVj = iVj_, where No- 
is the number of Fermi atoms with spin-cr, and that both Fermion components have the 
same single-particle energy. As the threshold energy of the Feshbach resonance is lowered, 
the character of the phase transition is found to continuously change from the BCS-type to 
BEC-type. In the latter regime, T c is strongly suppressed compared with the usual mean- 
field BCS theory and it approaches T c = 0.218Tp in the BEC limit. (Here Tf is the Fermi 
temperature of the non-interacting spin- up Fermi gas.) We show that the stable Bosons 
which appear as we enter the BEC regime are a strongly hybridized mixture of Feshbach 
molecules and preformed Cooper-pairs. The Feshbach component of these Bosons dominates 
the excitation spectrum in the BEC limit. 

In this paper, we extend our previous work |11]] by including the effect of a harmonic 
trap potential. This extension is clearly necessary since all experiments on Fermi gases are 
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done in some sort of trap potential. So far, the trapping potential has been examined within 
the mean-field BCS theory using the simple local density approximation (LDA). [IIJ In this 
paper, we also use the LDA but go past the mean-field approximation to include the effect 
of particle-particle fluctuations. We show how the BCS-BEC crossover can be observed 
from characteristic changes in the static atomic density profile, easily measured by standard 
techniques. 

Besides the phase transition temperature T c , it is also important to understand how 
strong-coupling fluctuation effects enter above T c . In particular, it is an interesting problem 
to clarify at what temperature the stable Feshbach molecules and preformed Cooper-pairs 
in the BEC regime first appear in the normal-fluid phase above T c . Our calculations show 
that above T c , there exists a characteristic temperature T* which describes the crossover 
from a free Fermi gas to a molecular Bose gas. 

In this paper, our treatment of the particle-particle fluctuations assumes that their cou- 
pling to the Feshbach resonance is weak and can be thus treated perturbatively. Most of 
our discussion is for this case, in which the BCS-BEC crossover is driven by decreasing 
the value of the Feshbach resonance threshold (2v) relative to twice the bare Fermi energy 
(2£p). When the Feshbach coupling strength is very large, as Kokkelmans and coworkers 



have found, |Tj| it is not clear that our simple treatment of the particle-particle fluctuations 



is valid. However, in this limit, one finds that the BCS-BEC crossover occurs for values 



2v ^> 2e-p because of the broad Feshbach resonance (see also Ref. [fTT|). Effectively, one is 



in the strong-coupling regime discussed by Nozieres and Schmitt-Rink, [|12| where the BEC 



phase is associated with preformed Cooper-pairs. When the threshold energy 2v becomes 
comparable to 2ef, our calculations again show that the Feshbach resonance takes over as 
the dominant component of the BEC phase. 

This paper is organized as follows: In Section II, we review the coupled Fermion-Boson 
model used in earlier work and discuss the LDA extension of the results we presented in Ref. 
[ n | . In Section III, we use this opportunity to expand on certain aspects which were not 



discussed in detail in Ref. |Tl| due to lack of space. The BCS-BEC crossover in a trapped 



system is discussed in Section IV. In Section V, we consider the character of the composite 
Bosons. The crossover from a free Fermi gas to a gas of long-lived molecular Bosons above 
T c is discussed in Section VI. Section VII considers the case of a strong Feshbach coupling. 



II. FORMULATION 



A. Coupled Fermion-Boson gas in a trap 



The gas of interacting Fermi atoms coupled to a Feshbach resonance can be described 
by the coupled Fermion-Boson model given by []7|,|^,|9| Jl^ , pTlJT8| , ^9|j20|| 



n = J2e p clc pa + J2(E° q + 2v)blb q 

per q 



p.q 



Here c pa is the creation operator of a Fermi atom with the kinetic energy e p = p 2 /2m. We 
assume that the system involves two atomic hyper-fine states which will lead to Cooper-pair 
Bosons and hence superfluidity. We describe these atomic states for simplicity using the 
pseudo-spin language a =],[■ The quasi-molecular Boson of momentum q associated with 
the Feshbach resonance is described by the operator b q . E® = q 2 /2M is the kinetic energy of 
this Boson and 2u represents the bottom of this resonance Bose spectrum, to be referred to 
as the threshold energy of the Feshbach resonance. The coupling to this Feshbach resonance 
Boson is described by last term in eq. ( |2.1| ) with a strength g T , in which two Fermi atoms 
form one 6-Boson and the 6-Boson breaks into two free Fermi atoms. The Hamiltonian also 
includes a standard BCS pairing interaction with the coupling constant —U < 0, originating 
from non-resonant interactions between two Fermi atoms (cr =1,1). Even assuming that 
this direct attractive pairing interaction U is weak, the effective pairing interaction strength 
can be strong due to the effect of the coupling to the Feshbach resonance. In this paper, 



we do not discuss the question of how good the model Hamiltonian in eq. (|2.1|) is in the 
context of ultracold Fermi gases, but simply refer to the literature []7|,|8|,|9| JlO| , p!6| . 
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Since in this model, one Feshbach-induced 6-molecule consists of two Fermi atoms, we 
have M = 2m and also impose the conservation of the total number of particles as 



pa 



= N F + 2N B . (2.2) 
As discussed in our previous paper, |TT[ this constraint is incorporated into the model Hamil- 



tonian in eq. fl2.1|) by replacing 7i — > 7i — fiN, where /i is the chemical potential. The re- 
sulting Hamiltonian has the same form as eq. fl2.1|) , but the kinetic energies of the Fermions 
and Bosons are now given by e p — ► e p — \i and E® — > E® — 2/x, respectively. 

Now we introduce the effect of a harmonic trap. In this paper, we consider an optical 
trap in which the two hyperfme Fermi atom states feel the same trap potential given by 

VWpO) = \ Yl m [ U lx x ) + ^lyy] + U 0z z j} , ( 2 - 3 ) 
Z 3 

where = (xj,yj,zj) is the position of the j-th atom. In our model, we also assume that 
the Feshbach molecules also feel this potential with m being replaced by the Boson mass 
M = 2m. We could relax this last assumption to deal more realistically with the dipole 
moment of the composite molecule but we leave this extension to future studies. Treating 



the resulting trapped Fermion-Boson gas within the LDA amounts to replacing the 

chemical potential \x by 

//(r) =/x-W r )- ( 2 - 4 ) 

The constraint on the total number of Fermi atoms is then replaced by iV = J dr[A^p(r) + 
2Nb(?)], where iVp(r) and A^b(i") are the number density of Fermi atoms and Bose molecules, 
respectively. 

The atomic hyperfme states may feel different trap potentials in a magnetic trap. This 
situation is similar to electron spins in superconducting metals under a magnetic field. In the 
theory of superconductivity, it is known that a magnetic field suppresses superconductivity 
and, under certain conditions, produces a non-uniform kind of superconducting state called 
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the Fulde-Ferrell state. p3| , |2^ . p5| The effect of fluctuations on this kind of Zeeman splitting 
in atomic Fermi gases is considered in another paper. |26|j 



B. Equations for T c : Thouless criterion 

The superfluid phase transition temperature T c can be conveniently determined using the 
Thouless criterion, in which the phase transition is characterized by a singularity in the 
particle-particle scattering vertex function r at uo = q = 0, associated with the formation 
of Cooper-pairs. [^7],[^| The four-point vertex function T within the generalized t-matrix 



approximation in terms of U and g r is given by the diagrams in Fig.l. In this figure, the first 
line describes ladder processes associated with the attractive interaction —U < 0, processes 



familiar in strong coupling superconductivity in electronic systems. |12| On the other hand, 
the second line includes the effective Fermion-Fermion interaction mediated by the 6-Boson 
propagator Dq defined in eq. ( |2.6| ), scattering processes which only arise in atomic Fermi 
gases coupled to a Feshbach resonance, such as in eq. Q2.1|). The equation of the vertex 
function T is given by 

r (p+, -p-;p' + , -v'-\ r) = -(U - g^D (q; r)) 

+ (U- g*D (q; r))i £ Go(p'|; v)G{-p"_ ] r)r(p'j, -p'L;p' + , -p'_; r), (2.5) 

P d" oj" 

where p + = (p + q/2,iu m + iv n ), p- = (p — q/2,iu; m ), q = (q,iu m ) and (3 = 1/T is the 
inverse of temperature; iu m and iv n are the Fermion and Boson Matsubara frequencies, 
respectively. Go and Dq represent the bare one-particle Fermion and Boson thermal Green 
functions given by 

G (p,iu m ;r) 



iuj m - e p + /i(r) ^ g ^ 



iv n ~ E° q -2v + 2fi(r) 

Since the modified chemical potential /i(r) includes the trap potential, the LDA Green 
functions and vertex function T also depend on position r. Multiplying both the sides of eq. 
by jjG (p + ; r)G (—p^; r) and summing over u m and p, we obtain 



1 n i " \nt " \tv n hi i \ (U — g^D (q;r))Il(q;r) 

a 1^ G o{P+, r )G{-p_;r)T(p + ,-p_;p -p_;r) = — — — — — — -. 2.7 

P St» 1-(U - g?D {q; r))II(g; r 

Here H(q; r) = II(q, iv n ;r) is the well-known particle-particle correlation function of the 
Cooper-pair-field operator defined by A(q, r;r) = J2 P c - p +«j/2|(t; i")c p + 9 /2t( r ; r ) in the ab- 
sence of U and g T : [ |Ill , |i1^ , [2l^ , j30f 



n(q,^ n ;r) = /* rfre^(T T {A(q, r; r) At(q, 0; r)}) 
Jo 



= \H G o(p + q/2, iw m + ii/,,; r)C (-p + q/2, -za; TO ; r) 

_ y 1 - /(e P +,/2 - Mr)) - /(gp-,/2 - Mr)) (2g) 
p + £ P - g /2 - 2/i(r) - ii/„ 

where /(e) is the Fermi distribution function. Substituting eq. (|2.7|) into eq. ( |2.5|) , we 

obtain 

l(q,z^,rj- 1 _ p7 _^ o(q)il/B;r)]n(qjil/B . r) - ^ 
We shall see that the effective attractive pairing interaction is given by 

U cS (q, iv n \ r) = U - g^D (q, iv n \ r). (2.10) 

According to the Thouless criterion, T c is given as the temperature at which the vertex 
function r(q, iv n \v) in eq. Q2.9| ) develops a pole at q = v n = 0. Within the LDA, since T 
depends on position r, this condition is satisfied at different temperatures T (r) depending 
on r. In this approximation, we have to regard the highest T (r) as the phase transition 
temperature of the trapped gas. |yj Since the density of particles is maximum at the center 
of the trap, we have T c = To(r = 0). Then the equation for T c within the LDA is given by 

l = f/ cff n(0,0;r = 0), (2.11) 

where 

U cS = U cS (q = 0,tu n = 0;r = 0) = U+ - 9 \ (2.12) 

2v — 2/i 

is the static part of the effective pairing interaction in the long wave length limit at r = 0. 



Equation (|2.11|) indicates that even if the direct pairing interaction U is weak, the effective 



pairing interaction U e g in eq. fl2.12|) can be large if g r is large or when the chemical potential 
of the 6-Boson (yU# = 2/i) approaches the bottom of the the Bose band 2v. More explicitly, 
we note that eq. ( |2.11| ) can be written in the well known BCS form 

1 = (U + g> 1 ) V tanh( ^-^ )/2Tc . (2.13) 
2v — 2fi p 2e p - 2\x 



As shown in the literature, p7|)P8| the pole which develops at T c at q = v n = signals 
an instability of the normal Fermi gas. At T < T c , the amplitude of this two-particle bound 
state can be shown to grow exponentially with time. This instability is removed by the 
appearance of a new phase below T c , described by the Cooper-pair order parameter. 

C. Fluctuation contributions to the chemical potential 

The chemical potential /i in eq. ( |2.13| ) is determined by the equation of the total number 
of Fermi atoms where a fe-Boson also counts as two Fermi atoms. In the weak-coupling BCS 
theory in the absence of the fe-Boson, this equation simply gives \l — e-p (where Ep is the 
Fermi energy in the absence of any 6-Bosons), assuming that the temperature dependence of 
/i can be neglected. Indeed, in most metallic (weak-coupling) superconductors, the deviation 
of n(T c ) from Ep is negligibly small as S/j, ~ {T c /e-f) 2 1. |32| On the other hand, when 



g T is large or the Feshbach resonance strongly enhances the effective pairing interaction 
Ucs, in analogy with the "strong-coupling" superconductivity discussed by Nozieres and 



Schmitt-Rink, |L2| the chemical potential \i is expected to deviate remarkably from e-p due 
to the effect of fluctuations associated with the Cooper-pair channel. We now discuss these 
fluctuations. 

The equation giving the number of Fermi atoms determines the chemical potential \i. 
This is obtained, as usual, |12| from first calculating the thermodynamic potential Q = 
/ drQ(r) and then using the identity iV = —dQ/dfi. The fluctuation contribution to the 
local thermodynamic potential fi(r) is shown in terms of diagrams in Fig. 2. In Fig. 2(a), 
we show the usual diagrams associated with the direct pairing interaction —U. Working to 
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all orders in U, this gives p"2|.|T3 



6Qu{r) = - £ e lUnS Ml ~ UU{q, zV n ; r)]. 
P 



(2.14) 



This is the contribution discussed in Ref. |12|]. In Fig. 2(b), we show the new fluctuation di- 
agrams contributing to fi( r ) associated with the Feshbach coupling interaction g r . Summing 
up to all orders, this gives 



«Ur) = ^E eW M 1 + 9r 2 A(q,^;r)n(q,^;r)] ) 



(2.15) 



where 



n(q, %v n \ r) = 



Il(q,iz/„; r) 



(2.16) 



1 - C/n(q,iv n ; r) 

describes the particle-particle bubble including the effect of U in the ladder approximation. 

We note that the total fluctuations contribution to fi(r) is the sum of the two contribu- 
tions in eqs. ( |2.14j ) and ( |2.15| ), which can be combined to give 



5Q(r) = - ]T e iv « 6 In [[1 - UU(q, iv n ; r)] [1 + g*D (q, iv n] r)fl(q, iv n - r)] 



J- e ^ 5 ln[l - C/^ (q, zV n ; r)n(q, zz/ n ; r)], 



(2.17) 



where the effective pairing interaction U e g(q, iv n \ r) is as defined in eq. ( |2.1U| ). In other 
words, the only effect of the diagrams in Fig. 2(b) is to renormalize the pairing interaction 
as U — > £/eff (q, iv n ') r ) i n ec L- ( ^-14| ). We also note that the collective fluctuations in eq. ( |2.17| ) 
(which are Bosons) are precisely those which are associated with the poles of the particle- 
particle vertex function T in eq. ( p.9|) . The thermodynamic potential for the contribution 
of the free Fermions and 6-Bosons is given by the sum of [B^] 



P 



n » = ■=Ec w lni3 - 1 (p I «/ B ;r). 

<^ ;»./'■. 



(2.18) 



Combining all these results and using iV = —dQ/dfi, we find the result 
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N 
where 



dr 



iV» + 2iV» - ~ £ e^"|-m[l - U eff (q, ix/ n ; r)n(q, ix/ n ; r)] 



(2.19) 



- 2 e^ 5 G (q, iu m ; r) = 2 £ /(e F - /i(r)), 



-~ e iVnS D (q,iu n ; r) = £n B (£ g ° + 2v - 2/i(r)). 



(2.20) 



Here tib(E) represents the Bose distribution function. Equations ( |2.11| ) and ( |2.19| ) provide 
us with a closed set of two coupled equations, from which we can obtain /i(T, N) and T C (N). 
In solving the coupled equations numerically, we have to take into account carefully the 
convergence factor e lUnS in the last term in eq. (|2.19|) . In the appendix, we explain how to 



sum up the Matsubara frequencies in this term. Above T c , we need only solve eq. ([2.1 9|) 
for fi(T,N). Our numerical results will be discussed in Sections IV and V. In the next 
section, however, we first discuss the physics involved in treating the fluctuations within the 
approximation developed in this section. This will allow us to better understand the nature 
of the Bosons involved in the crossover to a BEC regime. 



III. PHYSICS OF COUPLED 5-BOSONS AND COOPER PAIRS 



The contribution arising from the interaction in eq. ( |2.19| ) can be described in terms of 
Bose fluctuations given by the solutions of 

1 = U cS (q_, uo + i8] r)II(q, to + i5; r), (3.1) 

where we made the usual analytic continuation from the imaginary Bose Matsubara fre- 
quency to the real frequency axis as iu n — > u + i5. While eq. (|3.1|) appears relatively simple, 
it describes a complex situation involving the formation of preformed Cooper-pairs and their 
coupling to the bare 6-Bosons in eq. (|2.1|) due to the Feshbach coupling interaction g T . The 
preformed Cooper-pairs and the 6-molecules are coupled to each other, which leads to hy- 
bridization and only one kind of composite Bose molecule. It is useful to try to disentangle 
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these effect, although there seems no unique way of doing this in such a coupled excitation 
problem. 

Denoting the number density of atoms given by the last term in eq. ( |2.19| ) as 5Nfl{t), 
we have 

^FL(r) = \ £ e-^UDl + U eS (q)^] (_L_), (3.2) 

where for simplicity we denote II(q, iv n \ r) = II, D (q, iv n \ r) = D , U e s(q, %v n \ r) = U e s(q), 
and we have used the identity dD /dfi = — 2.D 2 . The effect of the particle-particle fluctua- 
tions described by Il(q, iv n ; r) naturally gives rise to a renormalized 6-Boson propagator D 
with a self-energy given by 

S(q, iv n ; r) = -# r 2 fl(q, iv n ; r), (3.3) 
where II has been defined in eq. ( 2.16| ). Thus we have 



b-\^iu n] r) = D l +g 2 v ti 

= iv n - [E° q +2v- ^ r 2 n(q, iv n - r)] + 2/x(r). (3.4) 



The self-energy £ in eq. ( |3.3| ) describes the effect of the decay of a 6-Boson into two Fermi 
atoms and their recombination to a 6-Boson, treating the non- resonant interaction —U 
within the ladder approximation. This self-energy can be also obtained by summing up the 
diagrams shown in the third line in Fig. 1 (shaded bubble in the figure). 

One may verify that the first term in the square brackets in eq. (|3.2|) renormalizes the 
number of 6-Bosons, giving 



77? 1 - U P ft(Q II 3 trf 



e 







i-M?)n Pfx i-u cS ( q )u 

= 2N b (t). (3.5) 
We recall that in our interacting Boson-Fermion model Hamiltonian, we have [see eq. 
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N = N F + 2N B , 



(3.6) 



where, within the LDA, 



N F = /dr£(ct ff (r)c p(T (r)>, 

p,(T 

N B = J dvJ2(bl(rMv)) 



(3.7) 



give the number of Fermions and Bosons, respectively. As discussed in many body texts, Ab 
is defined in terms of the renormalized single particle 6-Boson Green function -D(q, iu n ;r). 
In fact, one sees that A B = N B = j drN B (r), where iVe(r) is defined in eq. ( |3.5| ). This 
means that A F in eq. ( |3.7|) is given by 



A F = A° + 2A C , 



(3.8) 



where 2Ac is obtained from the second term in the square brackets in eq. ( |3.2|) as 

Equation ( p.8| ) shows that the expression for 2Nq in eq. ( |3.9| ) is the natural definition of the 
number of Cooper-pair Bosons which arise in the interacting Fermion gas with an attractive 
interaction given by U c s(q). 



The preceding analysis leads to the following expression for given in Ref. [TT 



N = A° + 2N B + 2N C , 



(3.10) 



where N B is the contribution from the 6-Bosons and Ac is the contribution from the Cooper- 
pairs, both being strongly renormalized through the effect of the Feshbach resonance coupling 
g r . The problem with the decomposition given in eq. ( |3. 1Q ) is that when g T ^ 0, one must 
remember that the 6-Bosons and the preformed Cooper-pairs are strongly hybridized. As 
a result, they give rise to a single Boson branch described by the pole of D in eq. 
namely 



u q = E +2v - g r Il(q, iu n -> u q + iS; r) 



-2/x(r) 



(3.11) 
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One can understand the basic physics by going back to eq. ( [2.171) and noting that 



i-c/ eff (g)n = [i-[/n][i+^ n] 

= [l-UU]+g 2 r U- 



(3.12) 



u - {E° q + 2v - 2/i(r)) ' 

This form shows clearly that the zeros of 1 — U e g(q)H correspond to a hybridized spectrum 
associated with the preformed Cooper-pair [1 — UU = 0] and the 6-Bosons [u — (E® + 2v — 
2/i(r)) = 0] by the coupling g 2 U. The zeros of eq. ( |3.12| ) are seen to be given by 



[1 - UU] [u - [El + 2v- 2/i(r)]] + g 2 U = 0. 
One sees that the same hybridization physics is given by the poles of D, i.e., 

n 



(3.13) 



u-(E» + 2v-2 f i(r))]+gt- 



0. 



(3-14) 



1 - UU 

All these versions emphasize that the renormalized Bosons given by D^ 1 = are the hy- 
bridized version of the bare preformed Cooper-pairs and the bare 6-Bosons. One can still ask 
what is the relative "weighting" of these two components in the renormalized Bose spectrum 
given by the pole of D, as given by the values of A^ B and N c in eq. (|3.10|) . In Ref. [)Tl 



we 



calculated both these contributions as a function of the threshold 2v in the case of a uniform 
Fermi gas. 

Finally, we note that one can write 

1 - U eS (q)U = dJd, 1 - UMD^U 



= D [D 1 
= D D M \ 



- 9?}n 



(3.15) 



Putting this into eq. ( 2.19|) , one finds that the term 2N^ is cancelled out by the Do factor 
in eq. (|3.15 ) and we are left with 



N = J dr[N^r) + 2N M (r)]. 



(3.16) 



Here the number density of composite molecular Bosons Nm{t) is defined as 
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2N M {r) = -Wz lUn& i;^D M \ (3.17) 



where D M is defined in eq. ( |3.15|) , 



D M \q,tu n ; r) = iv n - [E° q + 2u + (U/D - ^ r 2 )n] + 2//(r). (3.18) 



Note that II enters here, not IT as in eq. (|3.4j) . This composite molecule is a true hybrid 
of preformed Cooper-pairs and 6-Bosons and cannot be associated with either component. 
However, one can verify that the poles of Dyi are identical to those of D as defined in eq. 
(p|) , in that [see eqs. (|3~lp and (gl5|) 1 



;i - U^D- 1 = D' 1 - (£. - ^ r 2 )n = D^. (3.19) 

-^0 



IV. TRANSITION TEMPERATURE IN THE BCS-BEC CROSSOVER WITH A 

FESHBACH RESONANCE 

A. Renormalization for U and g T 

In this section, we consider the case of an isotopic harmonic trap {ujq x = u>o y = ujq z = c^o) 
for simplicity. It is convenient to take Ep as the unit of energy (for one Fermi species). As 
for the unit of length, we take the radius Rf of a non-interacting Fermi gas at T = defined 
by sf(Rf) = £f~ mu)\R\l1 = 0, or R F = 2eF/mujQ. We also take, as the unit of the number 
of particles, the total number of Fermi atoms given by 

Ar 4 [ RF 2 j f Mr) 2, (mUJofRF t A ,N 

N = — r 2 dr / jrdp = - — -, (4.1) 

7T Jo jo 24 



where Pf(?") = y2mep(r) is the Fermi momentum at r within the LDA. With this normal- 
ization, since iV and Rp always appear with the coupling constants U and g r as U = UN/ R F 
and gf = g^N/Rp, they can be absorbed into these coupling constants. (In the following, 
we simply write U and g v as U and g T .) In our numerical calculations, we simply intro- 

2/2 

duce a Gaussian cutoff e~ £ p /UJc , with u> c = 2e-p, in the momentum summation involved in 
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the particle-particle response function II(q, iv n ;r) defined in eq. fl2.8|) and thus also in eq. 

(EH)- 

The interactions —U and g T used in this paper should be regarded as renormalized quan- 
tities effectively involving the effects of screening by (pseudo-spin and density) fluctuations 
and renormalization in terms of high energy processes: Gor'kov et al. showed within the 
mean-field BCS theory that in the dilute limit, T c is suppressed by an effective interaction 
mediated by fluctuations described by particle-hole (p-h) bubble diagrams. |34| This effect 
was also studied recently by Combescot |35| in connection with the possibility of high T c 
superfluidity in 6 Li. This reduction of T c can be physically understood as a screening effect 
by (pseudo-) spin and density fluctuations, which weakens the pairing interaction between 
Fermi atoms. One should understand that —U and g x in our model in eq. (|2.1|) include this 
p-h screening. 

As for the renormalization of —U and g r by high energy processes, we note that our 
Fermi atomic gas has no physical cutoff-energy, in contrast to superconductivity in metals, 
where the lattice phonon Debye frequency acts the effective cutoff in the BCS gap equation. 
However, a prescription about the renormalization of the high energy interactions is implicit 
when we introduce a cutoff u c in the momentum summations. To see this, it is convenient 
to introduce an energy cutoff E c . This cutoff energy is assumed to be much higher than u c 
and T, and it may be arbitrarily large. We can then divide the summation in the equation 
for T c in eq. (|2.13|) into a low energy part e p = [0, u c ] and high energy part e p = [u c , E c ]. In 
this decomposition, when E c ^> u c > T c , we can take tanh(e p — /i)/2T c — > 1 in the 

high energy part. Since this replacement corresponds to ignoring the influence of the Fermi 
distribution [recall that in eq. (|2.13|) , tanhx/2T = 1 — 2f{x)\, the high energy region can 
be described as two atoms interacting in a vacuum. The effect of the surrounding gas and 
Fermi statistics is absent. Then the decomposition of eq. ( |2.13| ) is reduced to |36| 



i=^e tanh(£p : /i)/2rc +^ E I , (4-2) 

where U c s = U + g^/(2u — 2/x). This equation can be rewritten as the form 
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u, 



cff 



tanh(e p — /i)/2T c 



2e p — 2/i 



[0,w e 



tanh(e p — /x)/2T c 
2e p — 2/i 



(4.3) 



In this equation for T c , high energy two-particle scattering processes in the region [u c , E c ] has 

been absorbed into the renormalized coupling U^(u> c ) = U R (u c ) + g R2 (u c ) / (2v R (u c ) — 2/x) 

describing the low energy physics in the region [0, u c ]. The renormalized interaction U r (uj c ), 

Feshbach coupling g R (uj c ) and threshold energy 2u R (u c ) are thus given by 

U 



U R (u c ) 



2u R (uj r 



9? 

2 H[uc,Ec\ 2e p -2^i 

~-2v-g r 



(4.4) 



1-t/Er, 



>[uj c ,E c ] 2s p -2n 

Physically, U R (u) c ) is the non-resonant Fermion-Fermion interaction enhanced by high energy 
scattering processes and is diagrammatically described by Fig. 3 (a). The renormalized 
Feshbach coupling g R (u c ) involves the three point vertex correction coming from the high- 
energy processes in Fig. 3(b). Finally, the renormalized threshold energy 2u R (u c ) originates 
from the self-energy correction caused by the break of a Bose molecule into two Fermi atoms 
with high momenta, which can be described by the diagrams in Fig. 3(c). 

We should determine the renormalized quantities experimentally. In this regard, the 



renormalized quantities for the low energy limit (uj c = 0) can be introduced. |IE,l7| One can 
write the renormalized parameters in eq. (|4.4j) in terms of their low energy values (u c = 0) 
as follows: 



9?M 



2v R U) 



1 + U R (0) E[o,. c ] 

gK(0) 

1 + u R (0) e [0 ,. c] =V 

E 



(4.5) 



2u R (0) + g R (0) 



[0>wc] 2e p -2[i 



1 + ^ B (0)E [0) . C ]^ 

Written in terms of the renormalized variables given in eq. (|4.5|) , the BCS gap equation in 
eq. (|4.3|) no longer involves E c explicitly. Although we simply write the couplings as —U, g r 
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and the threshold energy as 2v in our model, we can regard them as renormalized quantities 
incorporating the high energy two-particle scattering processes, as discussed above. 

The present way of renormalization can reproduce the cutoff-free theories in Refs. 



I^ , p^5|j37| . When we write eq. ( |4.3[ ) using the renormalized variables in terms of the low 
energy value U^(uj c = 0) as given by eq. (|4.5|), we obtain 



t>5(0) En 



E 



'[O.Wc] 2cp — 2/i [0,w c ] 

This equation can be in turn rewritten in the form 



tanh(£ p — ji)/2T c 
2e p — 2/i 



(4.6) 



i = i/5(o) e [ tanl t p :" )/2r ° 



[0,w c ] 



2e: p — 2/i 2e p — 2/i 

rr* cm r tanh ( £ p ~ Z 1 )/ 2 ^ 1 



(4.7) 



[0,oo] " 2 ^ 2£ P ~~ 2 ^ 

Since the two terms in the square brackets almost cancel out each other in the high energy 
region [u c , oo] (we assume that uj c ^> T), we can eliminate the finite upper limit ui c in the 
summation. This cutoff- free expression is similar to the one used in Ref . . The cutoff-free 
BCS equation for T c used by Randeria et al. [[UjiEOl can ^ e a ^ so obtained from eq. fl2.13|) . 
Writing it as 



[0,E C ] 



tanh(e p — fi)/T c 



2e r , 



1 

2ei 



1 

2el 



(4.8) 



we can turn this into 



u cS E 



-tanh(e p — /i)/T c 1 



2e ri 



[0,E C ] 'I' I' 

where the renormalized interaction U e s = U + g^/(2i> — 2/i) is defined by 



(4.9) 



U 



i ' 



1 -UJ2[0,E C ] 2ip" 
1 

& = 9t 



1 -^E[0,E C ]2^' 



2v = 2v 



(4.10) 



^E r 



[0,E C ] 2^ 
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All these expressions for the equation for T c give physically the same result, assuming that 
one takes the appropriate values for renormalized parameters. 

With regard to the renormalization of the coupling constants, recent work by Kokkelmans 



et al. [01 has discussed the coupled Fermion-Boson model given in eq. ([2.1|) and incorporated 



the full two-body scattering theory into an improved fashion. Milstein et al. [T1J have more 
recently derived identical equations to eqs. (4) and (6) in Ref. [0 (these correspond to 
eqs. ( |2.13|) and ( [2.19| ) for a trapped Fermi gas), which include the effect of particle-particle 



fluctuations and the appearance of Cooper-pairs. They have extended our analysis in Ref. 
|TT| by renormalizing these equations in a way which treats the two-body scattering and 
bound states correctly. Their work thus leads to explicit numerical values for the various 



parameters in the renormalized model corresponding to eq. ( |2.1| ). 

In the rest of this Section and Sections V and VI, we discuss our results for T c and the 
region above, using U = 0.3ef and g r = 0.6£f- Our diagrammatic approximation for the 
fluctuations, as summarized in Figs. 1 and 2, assumes that the interactions are weak. In 
Section VII, we discuss the analogous results predicted by our equations ( |2.13| ) and fl2.19| ) 



using a value g T 3> £f, such as found in Refs. P,|I7|]. In this case, the Feshbach resonance is 
very broad and thus can have a huge effect on T c even if the threshold 2v ^> 2e f . 



B. Crossover behavior at T r within LDA 



Our numerical solution for T c and fi(T c , N) from eqs. ( |2.13| ) and ( [2.19| ) for a trapped 



Fermi gas gives results which are very similar to those for a uniform Fermi gas obtained in 



Ref. [Oj]. Figure 4 shows the BCS-BEC crossover in a trapped Fermi gas. As the threshold 
energy 2v is lowered, Fig. 4(a) shows that T c deviates from the mean-field BCS theory 
('BCS' in the figure) and approaches the BEC phase transition temperature of a gas of N/2 
free Bosons of mass M = 2m ('BEC in the figure). The chemical potential /i also changes 
remarkably in this crossover as shown in Fig. 4(b). As expected from the mean-field BCS 
theory, we obtain \i ~ ef for 2v > 2ef (BCS regime) while \x deviates from ef and approaches 
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v for v < (BEC regime). The latter result reflects that the 6-Bosons become the dominant 
contribution in the BEC regime. As a result, the phase transition occurs when the chemical 
potential of the 6-Boson /xb (= 2/x) reaches the bottom of the Bose band (2v) at the center 



of the trap (within our LDA). We note that we can rewrite the equation for T c in eq. ( |2.11| ) 
in the form 

2z/ = 2/j + # r 2 n(0,0,r=0). (4.11) 

This equation reduces to 2v = 2/x (= /xb) in the BEC limit, where the Fermi atoms and 
hence the particle-particle fluctuations described by II are absent. Since, in this BEC regime, 
the dominant excitations are N/2 free Bosons, we find the usual result for a trapped free 
Bose gas T c = u [N/2({3)} 1 / 3 , @ where uj is the trap frequency and ((x) is the usual zeta 
function (we set % = 1). Using the fact that T F = (QN^) 1 ^uj$ is the Fermi temperature for 
a trapped free Fermi gas (with iVj- = N/2), we can write this expression for T c in terms of 
T F to give 

T C =(^) 1/3 . F = 0.518T F . (4.12) 

Figure 4(a) shows that this T c is obtained in the BEC regime as the limiting maximum value 
in the trapped Fermi gas. We recall that Fig. 4 is for U = 0.3£ F and g r = 0.6e F . 

We note that this upper limit for T c is much higher than that of a uniform Fermi gas in 
the absence of a trap. JTT| In the uniform system, T c in the BEC regime is given by |TT|p!5 



T c = 2( — — ; N ) 2/3 £ F = 0.218T F , (4.13) 
where the Fermi temperature in the uniform system is given by T F = (67r 2 iV T ) 2 / 3 /2m. H 



One reason for this different maximum T c 's between the two cases comes from the different 
definition of T F . In addition, one should note that the harmonic trap potential enhances 
T c through the effective density of states (DOS) of Bose molecules. When we carry out the 
spatial integration in the equation for the total number of trapped particles in the BEC 



regime [N = 2 J driV B (r)], this equation is reduced to |J 
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poo I 

1 = 3 1 dEE-^— x (T = T C ), (4.14) 

where (3 = Tp/T and the energy in the integral is normalized in units of e-p. The corre- 
sponding equation in the uniform system is given by 

l = 3rf£'dEVE^- l . (4.15) 

This shows that the low energy region of the DOS has larger weight in the uniform gas 
(DOS oc VE) than in the trapped one (DOS oc E 2 ). As a result, the transition temperature 
(relative to Tp for each case) must be higher in the trapped gas than in the uniform gas. 



V. LONG-LIVED COMPOSITE MOLECULAR BOSONS 



A. Damping of 6-Boson and condition for stable molecules at T c 



The decay of 6-Boson into two Fermi atoms associated with the Feshbach resonance is 
described by the imaginary part of the (analytic-continued) self-energy of the renormalized 



Bose Green function in eq. (|3.4j) 



7 = — ImE(q, iu n — > u + i8;r) = g 2 lmll(q, uj + i5; r). (5.1) 

Since Jmll(q, u + iS; r) is given by 

mJm cosh f (% + J(u + 2/i(r) - E°)E°) 

ImU(q, u + iS; r) = ™ V ™ e(«; + 2^(r) - E°) In f 2 V , ^ = t q \ 

4tt/5,/^ qJ cosh f(f - J(u + 2(i(r) - E°)E^ 



2 \ 2 V v r ^ ' 1' 1> 

(5.2) 



where 0(x) is the step function, the damping rate 7 is finite only in the region u > uj^ = 
E® — 2/i(r). This step function restriction is clear when we note that the denominator of eq. 
can be written as [iu n + 2/x(r) — E®] — 2e p . At q = 0, eq. ( p.2| ) can be simplified to 



M(q = 0,cj + z5;r) = ^-^Q(u + 2/i(r))X + 2/x(r) tanh (5.3) 

47T v 4 
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Since the energy to dissociate one 6-Boson (with center of mass momentum q) into two 
Fermi atoms is given by [e p+q /2 — ^(r)] + [s- p + q /2 — M r )L "the threshold energy u> t h is found 
to correspond to the minimum energy to destroy one 6-Boson with q = 0. 

The damping rate 7 leads to finite lifetime of 6-Boson, or equivalently gives a finite width 
to the peak in the spectral density of the (renormalized) Bose excitation spectrum given by 

p B (q, r) = ImD(q, iv n — > 00 + iS; r). (5.4) 

7T 

In the BCS regime, where 7 is always finite for uo > (see the inset in Fig. 5(a)), the quasi- 
particle 6-Bose spectral density shows a broad peak, as shown in Fig. 5(a). This 6-Boson 
thus has a finite lifetime, which is given by the inverse of the peak width, due to decay into 
two Fermi atoms by the coupling to the Feshbach resonance. As one approaches the BEC 
regime, cJth > is realized (as shown in the insets in Figs. 5(b) and 5(c)). In this case, when 



the 6-Boson pole given by eq. ( |3.11 



u q = [E" - 2/i(r)] + \2v - ^II(q, uj q + 16; r)], (5.5) 

appears below u^, its lifetime is infinite. It thus appears as a ^-function peak in the 6-Boson 
excitation spectral density pe given by eq. ( |5.4j ). We show this case in Figs. 5(b) and 5(c). 
In Fig. 5(b), the excitation spectrum is still accompanied by a incoherent part in high 
energy region (u > Wth), in addition to the undamped bound state pole below uj^. On the 
other hand, the incoherent part is seen to be almost absent in the BEC regime in Fig. 5(c), 
evidence that a stable 6-molecule dominates the two-particle excitation spectrum. 
The condition for such a stable 6-Boson is simply that [|ll|] 

2z> = 2v - £ r 2 n(q, uj q + i5; r) < [£j - 2/i(r) > 0]. (5.6) 



This means that the renormalized threshold energy defined by 29 in eq. ( p.6| ) must be 
negative for the stable undamped 6-Bose molecules to appear. QTTJ In the region u = q = 0, 



this condition can be approximately written as 2z> ~ 2v — ^11(0, 0; r) < 0. Since 11(0, 0, r) 
can be shown to be positive and also decreases as r increases, the condition 2v < is first 
satisfied at the center of the trap. 
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We can divide the number of renormalized 6-Bosons in eq. ( |3.10|) into the contribution 
from stable long-lived part (N^~°) and a Feshbach resonating contribution with a finite 
lifetime (denoted by N^ >0 ). In order to evaluate the contribution from the stable part 
iVg -0 , it is convenient to introduce the spectral representation of the renormalized Bose 
Green function as [ I2" | 

1 roo I 

ZXq, iu n ;r) = / dz lraD{(\,iv n — > z + iS;r). (5.7) 

7T ./-oo iv n — z 

Using this expression, the summation in terms of the Boson Matsubara frequency in in 
eq. (|3.10|) can be carried out to give 



1 f f°° 

N B = V dr / dzn B (z)ImD(q, z + i8; r). (5.8) 

Using this, the contribution from the stable poles (with no imaginary part) is given by 

r poles r°° 

K = / dr J2 / dz5(z - [E q - 2//(r)] " - P r 2 n(q, z + 16; r)]) 

q J -°° 
poles 

= J dr^Z(q;r)n B (a; q ). (5.9) 



Here, Z(q; r) arises from what is called mass renormalization and is given by 

Z( q ;r)-' S l +J? %^£). (5.10) 

The sum in eq. ( |5.9|) is over the solutions u q satisfying eq. ( |5.5| ). 

Besides this Feshbach component of the two-particle bound state given by eq. ( |5.9| ). 
there is a component from Cooper-pair associated with iVc as defined in eq. (|3.9| ). One can 
rewrite eq. ( |3.9| ) in the form 



We next use the fact that 



on = (r^7n) a? (5 - 12) 
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and that for u> = u q [given by eq. fl5.5|) 1, one has the identities D 1 = —g^fl and U e g(q) = 
II -1 . With these relations, the contribution Nq~° of undamped Bosons associated with 



Cooper pairing can be reduced to the expression [1 I 



AT° = |/ ^E ^^ ^QiW, (5.13) 



where Z(q; r) has been defined in eq. ( |5 . 1 0|) . There is also a contribution from scattering 
states and hence we have Nq = Nq~° + Nq, where Nq represents the scattering contribution 
associated with damped particle-particle (or Cooper channel) fluctuations. 

Figure 4(c) shows the change of the character of the excitations in the BCS-BEC crossover 
in the trap. As v is lowered and approaches £p, the Feshbach resonance molecular state 
starts to play a role and the Feshbach damped molecules (N^ >0 ) appear. These Bosons 
are replaced by the stable molecules {N^~°) around v = 0.18ep for the parameters used in 
Fig. 4(c). At the same time, a contribution from preformed Cooper-pairs (Nq~°) appear. 
In the BEC regime, stable Feshbach molecules become dominant and N^~° approaches 0.5 
(= N/2). The contribution from preformed Cooper-pairs are seen to again disappear in the 
BEC limit v < 0. The increase in the total number of molecules is, however, slower than 



that found in the uniform system. Wl\ This is simply because the stable molecules initially 
only begin to appear at the center of the trap, rather than throughout the system as in the 
uniform g as case.. 

Figure 6 shows the number of composite Bosons Nyi = / drNu{r) as a function of v 
[where Nm(t) is defined in eq. ( |3.17| )1. Although each component of Boson (iVg -0 , N^ >0 , 



Nq and Nq) shows a singular behavior around v = 0.185f in Fig. 4(c), Nm (= ° + 
N£ >0 + N^ =0 + Nq) is found to increase smoothly as the threshold energy 2v is lowered, as 
expected. 

As we have discussed in detail in Section III, the strong hybridization induced by the 
Feshbach coupling interaction g T means that there is only a single two-particle Bose branch, 
whose energy is given by the solution of eq. (|5.5|) . Thus N^~° and iV^ -0 only represent the 
relative spectral weights of the two components (corresponding to long-lived 6-Bosons and 
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stable preformed Cooper-pairs) of this composite Boson. This decomposition is only meant 
as a way of understanding the physics of the BCS-BEC crossover which arises in the model 
given by eq. ( |2.1| ). This interpretation is, however, very useful to understand the difference 
between the present BCS-BEC crossover and that in strong-coupling superconductivity with 
no Feshbach resonance: In the latter case, the preformed Cooper-pairs are the dominant 
excitations in the strong coupling BEC regime. [|T^,|l4|,|]l| In the present case, for the values 
of the parameters U = 0.3£f and g r = O.Qef used in Fig. 4, although preformed Cooper-pairs 
make a contribution, the stable Feshbach molecules dominate the excitation spectrum in the 
BEC limit {y < 0). 

B. The binding energy of composite Bosons 

The binding energy of the long-lived composite Boson is the energy to break this molecule 
into two Fermi atoms. This dissociation energy can be evaluated from the energy difference 
between the bound state energy u; q in eq. (|5.5|) and the excitation spectrum of the free 
Fermi gas measured from the chemical potential /i. In particular, the binding energy of the 
zero-energy (u) q= o = 0) two-particle bound state with q = which appears at T c is given by 

^binding = 2[£ g=0 - /i(r = 0)] - u g=0 = 2\fl\ 

= 2|z?(q=0,o;, = 0;r = 0)|, (5.14) 

where we have used /i(r = 0) < and eq. ( |5.6| ) to obtain expressions on the right hand side. 
This result means that -Ebmdmg is equal to the absolute value of the energy of a renormalized 
6-molecule at the bottom of the renormalized Bose spectrum. 

The threshold energy co> t h(q = 0) of the two-particle spectrum is defined as where 
Jmll(q = 0, uj\ r = 0) given in eq. (|5.2|) becomes finite (see Fig. 5). We recall that ImU 
describes the excitation spectrum of the Cooper-pair channel fluctuations. This frequency 
u; t h(q = 0) is the same as the threshold energy of the continuum spectrum of the renormalized 
6-Boson self-energy in eq. ( p.3|) , since 7 = — JmE = g^lmX\[iv n — > uj + i5) oc ImU(uj + 18). 
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When a stable 5-Boson at u q=0 = is excited to a state with energy E > Uft, this excited 
6-Boson will decay into two Fermi atoms lading to a finite lifetime for this excited state. We 
conclude that u>th(q = 0) = -^binding, i-e., the particle-particle threshold energy o> t h(q = 0) is 
equal to the binding energy of a stable two-particle bound state. 

Since the chemical potential fj, approaches v in the BEC limit as shown in Fig. 4(b), 
we find from eq. (|5.14|) that ^binding — ► 2|i/| in this regime. This energy 2\v\ is just equal 
to the energy needed to transfer a free 6-Boson at the bottom of the Bose energy band 
(E® =0 + 2v = —2\u\) to two Fermi atoms with zero energy (2e q= o = 0). This result is 
consistent with the fact that the preformed Cooper-pair component in the bound state is 
absent in the BEC limit, leaving the 6-molecule component. 

The binding energy of the preformed Cooper-pair in the strong-coupling BEC regime has 



been also discussed in superconductivity, [|T^,|Tj.[n|.[l!| and it is useful to understand how 
the present results are related to this previous work. To see this, let us recall the cutoff-free 
equation for T c given by eq. Q4.9|). Since /i in the BEC limit, one can approximately 



take tanh(e p — /i)/2T = 1. Summing over p in eq. ( (4.9| ) , we obtain 



-(U + 9^—^r-) \ ■ (5.15) 



2u-2^ Air 

When we write the renormalized attractive interaction as U c ff = U + g^/(2u — 2/i) = 
— 47ra s /m, where a s is an s-wave scattering length, we find = l/2ma^, which is for- 
mally the same expression as that obtained in strong-coupling superconductivity discussions. 
|14] , ^5|J39| In the present model, where the pairing interaction mediated by the Feshbach res- 
onance is dominant in the BEC regime, we can neglect the non-resonant part U in U e g. 
Then in the case of v < 0, the solution of eq. (|5.15|) is f40 



128tt 2 |z/| 

1 + 



(5.16) 



v QAix 2 v m 3 gf 

The right hand side of this equation is reduced to unity in the BEC limit, v —>■ — oo, namely, 
the chemical potential for the 6-Boson (2/i) approaches the (renormalized) threshold energy 
2u, as one expects. 
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The relation between \i and the binding energy of a bound state in a two-particle system 
can be also obtained by extending the previous discussions given in the superconductivity 
case. 



I^JT5|J39[1 When we consider two Fermi atoms in a vacuum and employ relative coor- 



dinates, the scattering t-matrix for one Fermi atom with reduced mass m/2 and energy uj is 
given by 

t(uj) = (-u + g?^-) + (-u + g 2 T ^-) E -\-{-u + s r 2 ^r) + 

uj — lv' v uj — 2v' p uj — 2e p v uj — 2v' 

U + gh 1 



2u — U! 



(5.17) 



J P 2e p -uj 

The q = two-particle bound state is given by the pole of T(w), namely 

1=f/ e ff E^- ; (5-18) 

where U f eS = U + g% / (2u — uj) . When we renormalize the interaction in eq. (|5.18 ) as we have 
done to obtain eq. (|4.9|), eq. ( |5.18| ) is reduced to 

where U' cS = U + g^/(29 - co), in which the renormalized variables are given in eq. ( [4 . 1 0| ) . 
This equation is the same as eq. ( |5.17 ) in the BEC limit (/i 0, where tanh(e p — [x)/2T — > 1) 
when we replace uj — ► 2\i. Thus the energy of the bound state is related to the chemical 
potential by the relation uj = 2|/i|, which is consistent with our earlier result in eq. (|5.14j ). 



VI. STRONG COUPLING EFFECT ABOVE T c 
A. BCS-BEC crossover effect on the atomic density profile 



As expected, the change of the character of the particles from the Fermion to Boson 
in the BCS-BEC crossover shown in Fig. 4(c) strongly affects the density profile of atoms 
in the trap. This is shown in Fig. 7(a). In the BCS regime, the density profile is spread 
out, as shown by result for v = £p. In this regime, the density profile is mainly composed 
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of Fermi atoms, and the Pauli exclusion principle between Fermi atoms effectively acts as 
a repulsive interaction. This pushes the atoms away from the center of the trap. As the 
system approaches the BEC regime {y < 0), the number of stable Bosons gradually increases 
and thus the repulsive effect of the Pauli exclusion principle becomes less important. As 
expected for Bosons, the particles are now seen to cluster near the center of the trap as v 
is lowered, as shown by Fig. 7(a) [see also the inset in Fig. 7(a) for v — 0]. These results 
indicate that the easily measured atomic density profile may be a very useful signature in 
looking for the BCS-BEC crossover phenomenon experimentally, especially as one gets close 
to T c . 

We note that the atomic density profile shows a steep decrease near the trap center at 
T c (Fig. 7(a)), but this is absent at T = 1.5T C (Fig. 7(b)). The cusp-like behavior at T c 
originates from the Boson component (N-q and Nq) in the density, as shown in the inset in 
Fig. 7(a). In the extreme BEC limit, since the gas is described by N/2 free Bose gas in a 
trap, the density of atoms is proportional to g 3 / 2 (z(r)), |38| where 

9»{z) = En=i z n /n 3 / 2 and 



the fugacity is z(v) = e^~ 2v ~ 2VtIa:p ^'' T . At T c , since the BEC is realized when the chemical 
potential of Bosons (2/z) reaches the bottom of the Boson excitation spectrum (2z/), we find 
z(r) — > 1 at the trap center. From the well known behavior of gs/2{z) at z = 1, the atomic 
density profile shows a finite slope at r = at T c . In contrast, this sharp cusp at the trap 
center is absent for temperatures above T c , where 2fi < 2u and hence the fugacity z(r) is no 
longer close to unity. We also show in the inset in Fig. 7(b) the density profile above 
T c when the temperature is fixed as T = 0.75Tp. In this case, the density profile becomes 
more spatially diffuse as v is increased. This is because of the increasing dominance of the 
Fermions and the resulting enhancement of the effect of the Pauli exclusion principle. 

B. Crossover from a Fermi gas to a Bose gas above T c 

It is useful to clarify when and how the stable Bosons {N^~° and Nq~°) at T c shown in Fig. 
4(c) disappear as the temperature is increased above T c . Figure 8(a) shows the temperature 
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dependence of the number of particles above T c (for v — 0). The number of stable Bosons 
(N^ =0 and N^ =0 ), as well as the finite lifetime Bosons (N^ >0 ) and the scattering contribution 
(Nq), are found to all decrease at higher temperatures. This result is also obtained when 
v < 0. Thus only the free Fermi atoms contribute at higher temperatures well above T c . 

The reason that the Fermi atom contribution dominates at high temperatures has a 
very simple explanation in the difference of the chemical potentials between a Fermi atom 
(/iF = /-0 and a 6-Boson (/iB = 2/z). As the temperature is increased, the chemical potential 
/i must decrease so as to conserve the total number of atoms [N— (Fermi atoms) +2 x (Bose 
molecules)], as shown in Fig. 8(b). Thus for T ^> T c , the Boson excitation energy E^+2u—2fi 
always increases faster than the excitation energy of the Fermions e p — //, simply due to the 
factor two in front of /i in the Bose excitation energy. As a result, the occupation of the 
Bose band always approaches zero for T ^> T c , even though the Bosons can be the dominant 
excitation at T c in the BEC regime. 

However, Fig. 8 still shows that there is a significant region above T c where the number 
of long-lived Bosons is substantial. It is thus useful to define a "crossover temperature" T* 
which separates the Fermi gas and the Bose gas regimes, even though the increase in the 
number of stable Bosons occurs continuously. We define T* as the temperature at which 
iVg -0 = 0.05, or when 10% of the Fermi atoms combine to form stable 6-molecules (Fig. 9). 
Figure 9 shows that T* rapidly increases when v is small or negative, for which case stable 
Bosons appear at T c (v/e F < 0.18, for the parameters used). Even for v ~ 0, T* ~ e F is 
considerably larger than T c . Since this "Bose gas regime" lying between T c and T* is fairly 
wide, one might expect to observe the formation of the stable molecules somewhat above T c 
if we can decrease the threshold energy 2z/ enough, namely to values where the superfluid 
transition at T c would be of the BEC type (see Fig. 4). 

Figure 10 compares the temperature dependence of the number of particles and the chem- 
ical potential (lines in the figure) for v = — e F with the case of a free (U = g T = 0) Fermion- 
Boson mixture (filled circles and squares). In the non-interacting Fermion-Boson mixture, 
we only impose the constraint that the total number of atoms is conserved [N= (Fermi 
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atoms)+2x (Bose molecules)]. Figure 10 shows that deep in the BEC regime (u = —£p), 
this constraint alone determines the temperature dependence of the number of particles and 
the chemical potential. As far as these quantities are concerned, the effect of including the 
interaction g T and U is clearly negligible. 



VII. THE CASE OF STRONG FESHBACH COUPLING 



In all of the numerical results presented so far, we have considered the case of a weak 
Feshbach coupling described by g r = 0.6ep. In this regard, recent calculations have argued 
that strong coupling (g r ^> e F ) was a more realistic description of Feshbach resonances 
observed in ultracold gases. P, |TE| , |r7|] In this section, we investigate what arises in this strong 
Feshbach coupling regime, although the present treatment of fluctuations implicitly assumes 
a weak Feshbach coupling. We present some numerical results for T c in the limit of a large 
Feshbach coupling (broad resonance) for a uniform gas. In this regard, we note that our 



analysis in the present paper and in Ref. [ 1 1] has shown that the BCS-BEC crossover in a 



trap is very similar to that in a uniform gas, at least for the case of weak Feshbach coupling. 

Figure 11 show the BCS-BEC crossover in the case of g r = 20ep ^> Sp. Figures 11(a) 
and 11(b) show that the crossover behavior already occurs around v = 150^ which is much 



higher than Ep. Results similar to this are presented in Ref. [17]. Since the bottom of the 
Bose excitation spectrum 2v is then still much higher than the Fermion band, the Feshbach 
related 6-Bosons described by N B = N^ =0 + N^ >0 are almost absent in this crossover regime, 
as shown in Fig. 11(c). The high-energy 6-Bosons only contribute to virtual scattering 
processes involving Fermions, which mediate the pairing interaction between Fermi atoms. 
Thus, when the Feshbach coupling g r is very strong, the character of the BCS-BEC crossover 
becomes identical to that in the case of strong-coupling superconductivity as discussed by 



Nozieres and Schmitt-Rink. [12| The fact that the crossover to the BEC region occurs at such 
large values of 2v is easy to understand. Since the chemical potential \i is always of order 
of £f or less, at large values of 2v and g T , we have U e s = U + g^/(2u — 2/i) ~ g^/2v. Since 
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the strong-coupling effect should be most pronounced when U e g becomes comparable to ef, 
the BCS-BEC crossover will occur at U e g ~ £p, which gives 2u ~ g x / £ F- For g T = 20ep, this 
predicts that the crossover will occur at v ~ 200£f, consistent with the numerical results in 
Fig. 11(a). 

That the case of a very broad Feshbach resonance (g r ^> e-p) is equivalent to the case 



studied in Ref. [|12l , |15l is also shown clearly by the fact that stable preformed Cooper-pairs 
are the dominant excitation [see Fig. 11(c)] in the crossover region at v ~ 150£f- The 
Feshbach 6-molecule component only becomes dominant when v < 0. This result is in 
contrast to our results for a weak Feshbach coupling g x = 0.6sf as shown in Fig. 4. In that 
case, the 6-molecule becomes the dominant component at the BCS-BEC crossover. 

Figure 12 shows the total number of composite Bosons N M = N^~° + N^ >0 + N^~° + 
as a function of u, for g r = 20ep- Although the individual Cooper-pair components Ng~° and 
Nq show [see Fig. 11(c)] singular behavior in the BCS-BEC crossover regime (at v ~ 150£f), 
A^m itself increases smoothly as the threshold energy is decreased. 

Since the BCS-BEC crossover around v = 150ep shown in Fig. 11 is intrinsically the same 
phenomenon as discussed by Nozieres and Schmitt-Rink, \T2\ the origin of the small peak in 
T c shown in Fig. 11(a) is the analogue of the one formed also in the case of strong-coupling 
superconductivity. [|T^,[l4|,|l5[] However, Haussmann |^l| showed that the slight peak in T c 
obtained at the BCS-BEC crossover in strong-coupling superconductivity disappears when 
one includes higher order scattering processes beyond the simple t-matrix approximation, 
such as used in the present paper. Thus, we regard the small peak T c in Fig. 11(a) as 
probably an artifact of our simple t-matrix approximation. As in the weak Feshbach coupling 
case shown in Fig. 4, we expect that the maximum transition temperature T c will again be 
given by the BEC expressions in eqs. (|4.12| ) and ( f4.13| ). 
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VIII. SUMMARY 



To summarize, we have extended our recent work [ I I ] and investigated the effect of 



particle-particle fluctuations in conjunction with an atomic Feshbach resonance on the 
superfluid phase transition in a gas of Fermi atoms in a harmonic trap, using a simple 
LDA approach. The BCS-BEC crossover predicts a maximum transition temperature of 



T c = 0.518Tp, somewhat larger than T c = 0.218Tp in a uniform gas. |TTJ In the BEC regime, 
the excitation spectrum is dominated by the stable long-lived Feshbach-related molecules 
in the case of a small Feshbach coupling parameter g r . This is in contrast with results for 



the BCS-BEC crossover with no Feshbach resonance, [pT2] , p~5|| where the preformed Cooper- 
pairs are the dominant bound states. For a very broad Feshbach resonance (g r ^> £p) as 
studied in Refs. P Jl6|Ji7|| , we find that the crossover region is also dominated by preformed 
Cooper-pairs. 

The BCS-BEC crossover is shown to strongly affect the atomic density profile in the trap, 
which can be measured. In the BCS regime, where Fermi atoms are dominant excitations, 
the spatial distribution of the atoms is spread out from the trap center due to the Pauli 
exclusion principle. The atoms increasingly cluster at the center of the trap as the threshold 
energy of fe-Bose excitation spectrum is lowered and the the number of stable composite 
Bosons increases. 

In the BEC regime, we show that stable molecules can exist above T c . As the temper- 
ature increases, however, the number of these Bosons decreases. We introduce a crossover 
temperature T* defined as the temperature above T c at which 10% of the Fermi atoms have 
formed long-lived 6-molecules. Plotting T* as a function of the threshold 2z/ gives a sort of 
"phase diagram" in the normal phase which may be useful in experimental searches for the 
BCS-BEC crossover. 

The experimental observation of Cooper-pairs in trapped Fermi gases will be a very 
exciting milestone in physics. What makes this topic even more interesting is the possibility 
that we can observe the BCS-BEC crossover in such systems. As we have shown in this 
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paper, this crossover involves very interesting physics and, moreover, it leads to characteristic 
changes in the properties of the trapped gas. The changes in the atomic density profile may 
provide a crucial experimental signature for the existence of Cooper pairing and superfluidity. 

We will discuss the BCS-BEC crossover in the superfluid region below T c in a future 
paper. In particular, we will discuss how the collective modes vary as we pass through the 
crossover region. |42f 
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APPENDIX A: CALCULATION OF FLUCTUATION CONTRIBUTION TO THE 

NUMBER OF FERMI ATOMS 



The density of atoms associated with fluctuations given by SN-p^r) in eq. ( 2.19 ) can be 



split into two contributions [using eq. ( 2.171) 1, |4 



5Ni 



fl r 



P 



(Al) 



Since the second term involving ^- ln[l + g 2 D fl] goes as v~ 2 in the large v n limit, the 
convergence factor e lu,lS has no effect. However, this convergence factor cannot be dropped 
in the first term in eq. ( |A1| ), because the quantity 1 _ U un ^ only goes as v~ l for large v n . In 
numerical calculations, it is difficult to sum up the Matsubara frequencies directly taking 
into account the convergence factor involving a infinitesimally small number 6. One method 
to avoid this difficulty is to employ a spectral representation as carried out by Nozieres and 
Schmitt-Rink. [|T^] However, this method requires high numerical accuracy in executing the 
resulting energy integration. We used a different approach. 

To calculate the first term in eq. (|A1|) [= 5Np^(r)] without employing a spectral repre- 
sentation, we compare 6N$(t) to the expression without the convergence factor, 



= - V- u 



(A2) 



This quantity is easy to calculate numerically. When we transform the summation over the 
discrete Matsubara frequencies into an integral in the complex plane in the usual way, we 
obtain 



5N, 



(i), 

FL I 



1 



(2), 
FL I 



2m 



2-ni 



dze Sz n B (z)- — — 

Co 1 — U\x{tu n 



dH(iv n 



E 



dzn B (z) — 
Co 1 



U 



z) d\i 
dU(iu n — > z) 



(A3) 



UYi{iu n — > z) d\i 

where the path Cq is shown in Fig. 13. we are allowed to add the circular contours C\ and 
C*2 for the integration in <WpL(r) without changing the result, because the contributions 
from C\ and C 2 vanish in the large radius limit [R — > oo) due to, respectively, the factors 
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ub(z) and e Sz . The integration in 5N^(r) is thus reduced to summing up the poles inside 
the closed contours Co + C\ and Co + C2 shown in Fig. 13. 

We next show how 5N^(t) is related to 5Np^(r) defined in eq. (|/^2|). First, we note that 

(2) 

when dealing with <5iVp L (r), the integration along the path C% gives a finite contribution 
when we add the paths C\ and C2 due to the absence of the convergence factor. However, 
if we write the integration in dNp^r) as 



dz = dz- dz, (A4) 

Co JCo+d+d J C'2 

we find that the first term can be calculated by summing up the contribution from the poles 
inside the closed paths Cq + C\ and Co + C2, and the result is identical to 5Np^(r). The 
last term in eq. (~K4) [= 5Np™ T (r)] can be evaluated to give, in the limit R — ► 00, 

SN™(t) = -— [ ^r ^(W2-MW) a/(£ p - g /2-Mr)) i 

= £/ V e- £ p^ gfetgZgJli^jl , (A5) 

Here, e~ £ p^^ is the Gaussian cutoff introduced in Section IV, and we have used eq. ( |2.8| ) 
and changed the integration variable using z = Re %e . Thus we have shown that SN^{r) 
involving the convergence factor is related to iVp L (r) without the convergence factor by 

6N®{r) = flvgg(r) - 6N^(v). (A6) 

Using this expression, we can calculate SN^^r) numerically by evaluating the simpler ex- 
pression 5NpQ(r). 

The essence of this prescription may be easily understood from the following example: 
Consider the two quantities: 



^ 2) (q) = -E ' 



%Vn E q (A7) 



Un iv n - E q 

In this case, the frequency sum in N#\q) can be easily done analytically and reduced to 
the Bose distribution function n-B[E q ). The sum in gives 
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„ ^ + £„ 2 2 2 



The contribution from C2 in the last term in eq. (|A4j) is given 

1 /■ 11 

JAr corr = / dz 



" 2iri Jc 2 :R^<x> z — E q 2 



Thus we can easily verify that = — 5N^™ , as in eq. 
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FIGURES 

FIG. 1. Generalized i-matrix approximation for the particle-particle scattering vertex function 
r(p + , —p^,p' + , —p'_). The solid and dashed lines represent the one-particle thermal Fermion Green 
function Go and Boson Green function Dq, respectively. The first line includes the ladder processes 
associated with the direct pairing interaction —U. The second line includes the effect of the 
Feshbach resonance with the coupling g T . The shaded bubble in the second line includes the ladder 
diagrams shown in the third line. 

FIG. 2. Fluctuation contributions to the thermodynamic potential fi. (a) Ladder diagrams 
associated with the direct pairing interaction —U. (b) Diagrams arising in the presence of the 
Feshbach resonance. The shaded bubble is the same as in Fig. 1. 

FIG. 3. (a) Renormalization of the non-resonant pairing interaction —U, as given in eq. (4.4). 
The solid line with solid circle means the Fermion Green function in the high-energy region [u c , E c ]. 
(b) Renormalization of the Feshbach resonance coupling g r , which is described by the three point 
vertex correction, (c) Renormalization of the threshold energy 2v associated with the self-energy 
correction. 

FIG. 4. BCS-BEC crossover in a trapped Fermi gas, for the parameters U/ef = 0.3 and 
9t/£f = 0.6. (a) The superfluid phase transition temperature T c as a function of v. In this figure, 
BCS is the result neglecting particle-particle fluctuations while BEC is the value for T c for a gas 
of N/2 Bosons of mass M = 2m. (b) The chemical potential [i at T c as a function of v. The 
dashed line (2/x = 2v) shows the chemical potential for a trapped gas of N/2 6-Bosons. (c) Change 
in the character of particles through the BCS-BEC crossover region. We note that scattering 
contribution Nq becomes negative below v/e-p — 0.18, where the stable Feshbach molecule (iVg - °) 
and preformed Cooper-pair (Nq~°) components first appear. The sum N^,~° + Nq is positive. 

FIG. 5. The spectral density pb of the renormalized Bose Green function D, for T = T c and 
r = and q = pp, where pf is the Fermi momentum for a free gas of N/2 spin up Fermions. In 
each panel, the inset shows the frequency dependence of the damping 7(0;). 
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FIG. 6. The number of Bosons Am [see eqs. (3.16) and (3.17)], as a function of the threshold 
energy v. The number of Fermi atoms A^p is also shown. 

FIG. 7. Atomic density profile at T c as a function of position in the trap in the BCS-BEC 
crossover, (a) T = T c . (b) T = 1.5T C . Inset in Fig. 7(a): Density profiles of iVg(r), iVg (r) 
and A^ =0 (r) at T = T c in the case of v = 0. Inset in Fig. 7(b): Density profile above T c when 
the temperature is fixed as T = 0.75Tf. These results are for U/ef = 0.3 and q t /ef = 0.6. The 
number density of atoms N(r) is related to the total number of atoms N as N = j dxN(r). One 
Boson counts as two Fermi atoms in this figure. We note that different temperatures are used for 
the three cases {v/ef = ±1,0) in Fig. 7(b). Since T c increases as v is lowered [see Fig. 4(a)], 
the temperature for v = ef in Fig. 7(b) is the lowest among the three cases plotted. As a result, 
the spatial diffusion of atoms in the density profile originating from the increase of the (averaged) 
kinetic energy of atoms due to a finite temperature is weakest in the case of v = ef and hence the 
density profile at the trap center is largest. 

FIG. 8. (a) The number of particles above T c when the threshold energy is zero, for U/ef = 0.3 
and g r /EF = 0.6. (b) Temperature dependence of the chemical potential fx above T c . 

FIG. 9. Schematic phase diagram for Bose molecule phase. For a given value of u, the temper- 
ature T*{u) is defined where the number of stable Feshbach molecules Ag _0 equals 0.05 AT. These 
results arc for U/ef = 0.3 and g r /£F = 0.6. Above T*(u), the system can be regarded as a gas of 
Fermi atoms, while below this temperature, one has an increasing fraction of pairing into stable 
molecules. The number of stable Feshbach molecules and the preformed Cooper-pairs at T c are 
also shown. 

FIG. 10. (a) The number of Fermi atoms and the stable Feshbach molecule component 
N^~° for temperatures above T c , for v = — ep. (b) The chemical potential \i above T c . The lines 
show the results for U/ef = 0.3 and g T /EF = 0.6. The circles and squares are the results for a 
non-interacting Boson-Fermion gas (U = g T = 0) but with the constraint that the total number of 
the particles A^ = [Fermi atoms+2xBose molecules] is fixed. 
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FIG. 11. BCS-BEC crossover in a trapped Fermi gas in the case of a strong Feshbach coupling 
<?r/ £ F = 20. We take U/ep = 0.3. (a) T c as a function of v. (b) The chemical potential \i at T c . 
(c) Change of the character of particles through the BCS-BEC crossover. The number of stable 
Cooper-pairs A^ -0 exceeds N/2 (= 0.5N) around the sharp peak shown in panel (c). However this 
is canceled by the negative value of the scattering contribution Nq. The total number of Bosons 
is always smaller than N/2 (see the graph of Nm in Fig. 12). 

FIG. 12. Change of the number of Bosons A^m in the BCS-BEC crossover in the case of strong 
Feshbach coupling, taking g T /ep = 20 and U/e-p = 0.3. In this figure, the contribution from N^ >0 
is negligible. Compare with the analogous results for weak Feshbach coupling in Fig. 6. 

FIG. 13. Path for the complex integration discussed in the Appendix. The dots show the 
imaginary discrete Bose Matsubara frequencies iu n . 
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